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ABSTRACT 

1ES 1959+650 is one of the most remarkable high-peaked BL Lacertae objects (HBL). In 2002 it 
exhibited a TeV 7-ray flare without a similar brightening of the synchrotron component at lower 
energies. This orphan TeV flare remained a mystery. We present the results of a multifrequency 
campaign, triggered by the INTEGRAL IBIS detection of 1ES 1959+650. Our data range from 
the optical to hard X-ray energies, thus covering the synchrotron and inverse-Compton components 
simultaneously. We observed the source with INTEGRAL, Swift/XKT, and UV-Optical Telescope, 
and nearly simultaneously with ground-based optical telescope. The steep spectral component at X- 
ray energies is most likely due to synchrotron emission, while at soft 7-ray energies the hard spectral 
index may be interpreted as the onset of the high-energy component of the blazar spectral energy 
distribution (SED). This is the first clear measurement of a concave X-ray - soft 7-ray spectrum for 
an HBL. The SED can be well modeled with a leptonic synchrotron- self-Compton model. When the 
SED is fitted this model requires a very hard electron spectral index of q ~ 1.85, possibly indicating 
the relevance of second-order Fermi acceleration. 

Subject headings: BL Lacertae objects: individual (1ES 1959+650) - galaxies: jets - radiation mech- 
anisms: non-thermal - X-rays: galaxies 



1. INTRODUCTION 

Blazars are very bright from radi o frequencies to 
gamma-ray energies (|Ulrich et allll997l ) due to the align- 
ment of their jets with respe c t to the line of sight o f 
the observer (|Antonuccil 119931 lUrrv fc Padov ani 1995). 
Relativistic beaming in these outflows ([Reesl I1966T) en- 
hances the observed flux and shortens the variability 
timescales. The spectral energy distribution (SED) typ- 
ically shows a non-thermal two-component structure. 
The low-energy emission is generally believed to be pro- 
duced through synchrotron radiation of relativistic elec- 
trons. In contrast, the nature of the high-energy com- 
ponent is still being debated. A possible explanation 
is that the same electron population that is responsible 
for the low-energy component generates the high-energy 
component through Compton sc attering (jSikora et al.l 
119941 : iDermer fe Sch lickciscr 1993J. Alternative expla na- 
tions are proton-initiated cascades (|Mannhefrnl I1993T ) or 
proto n-synchrotron emi ssion from ultrarelativistic pro- 
tons (jMucke et al.ll2003l ). 

The position of the synchrotron peak defines two 
classes of BL Lac objects: high frequency peaked BL 
Lac objects (HBL, peak at UV-X-ray frequencies) and 
Low frequency peaked BL Lac objects (LBL, peak in 
the IR-optical band) 1ES 1959+650, at z = 0.046 
(|Veron-Cettv fe Veronl 120061 ), is one of the best-studied 
member of the former group. It was first det ected at X- 
rays d uring the Slew Survey of Einstein-lPC (jElvis et al.l 
1992). Further observations detected the object at 
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X-ray s with ROSAT and BeppoSAX (|Beckmann et all 
2002). Based on the X-ray/radio versus X-ray/optical 
color-col or diagram the source was classified as BL Lac 
object bv lSchachter et a l. ( 1993) . The blazar was a lso de- 
tected at 7-rays by EGRET ([Hartman et al1ll999f ). Very 
recently, the source was included in the Fermi first cat - 
alog of active galactic nucleus (AGN fAbdo et aD r2010'). 
1ES 1959+650 was observed several times at TeV ener- 
gies (jAharonian et all 120031) . On 2002 June 4, followed 
by the detection of a strong TeV 7-ray flare of the source 
with the 10m Whipple Cerenkov Telescope, Target of 
Opportunity observations in optical and X-rays were per- 
formed. Despite an increased activit y in the 7-rays, no 
flux v ariation was detected in X-rays (|Krawczvnski et all 
120041 ) during simultaneous observations by the Rossi X- 
Ray Timing Explorer (RXTE). As a possible explana- 
tion, a ha dronic syn c hrotr on-mirror model has been sug- 
gested bv iBottcherl (|2005f ). bu t this suffers fr om rather 
extreme energy requirements (|B5ttcherl 120071 ). Hence, 
the origin of the orphan TeV flare is at present not 
understood. It is therefore important to monitor the 
synchrotron and inverse-Compton components simulta- 
neously. 

In this Letter, we report on a multiwavelength cam- 
paign on 1ES 1959+650 ranging from the optical to hard 
X-ray energies. In Section 2, we describe the observa- 
tions of the multiwavelength campaign, and in Section 3 
its results. In Section 4, we summarize our main conclu- 
sions. 

2. OBSERVATIONS 

2.1. The 2007 multiwavelength campaign 

1ES 1959+650 was monitored by INTEGRAL (from 
2007 November 27 th to 2007 December 1 st 2007 ) dur- 
ing the North Ecliptic Pole Key Program (proposal id: 
0531000). The source was detected by IBIS/ISGRI 
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TABLE 1 

Log of INTEGRAL /IBIS, Swift/VVOT and Palomar 60- inch Observations. 



Instrument 


Start Date UTC 


Exp (s) 


Band 


mag 


Err 


Flux 6 


UVOT 


2007-Nov-16 02:34:55 


1069 


u 


14.96 


0.02 


- 


UVOT 


2007-Nov-23 03:00:46 


479 


uvwl 


14.98 


0.03 


- 


UVOT 


2007-Nov-23 04:34:45 


539 


uvwl 


15.09 


0.02 


- 


UVOT 


2007-Nov-30 17:54:42 


98 


uvw2 


15.33 


0.03 


- 


UVOT 


2007-Nov-30 17:56:10 


165 


uvwl 


15.19 


0.03 


- 


UVOT 


2007-Nov-30 17:58:55 


388 


u 


14.78 


0.02 




UVOT 


2007-Nov-30 18:05:23 


318 


b 


14.63 


0.02 




60 


2007-Dec-03 01:56:02 


30 


9 


15.28 


0.03 




60 


2007-Dec-04 01:56:33 


30 


9 


15.29 


0.03 




60 


2007-Dec-05 01:54:03 


30 


9 


15.32 


0.03 




60 


2007-Dec-ll 02:25:35 


30 


9 


15.38 


0.03 




ISGRI 


2007-Nov-24 16:51:47 (624) 


55000 








4 -°2.6 


ISGRI 


2007-Nov-25 19:47:58 (625) 


198000 








3 8 1 ' 3 


ISGRI 


2007-Nov-28 20:16:36 (626) 


198000 








5 3 1 - 3 


ISGRI 


2007-Dcc-Ol 19:19:58 (627) 


208000 








2 fi 1 - 5 
z.o- lA 



"For^ISGRI observation, revolution number is reported in ''fluxes that arc computed in the 20 - 40 keV energy range in units of 10 



erg 



(|Bottacini et al.l 120071 ) in the 20 - 40 keV band in an 
active state at a significance of 7.2 a. Simultaneous X- 
ray and UV observations were obtained as a Target of 
Opportunity with Swift on 2007 November 30 th . From 
December 3-11, the Palomar (60) telescope performed 
optical photometry. Details of the multiwavelength cam- 
paign can be found in Table [TJ 

2.1.1. UV to Optical Observations 

UVOT. 1ES 1950+650 was observed in the u, b, uvwl, 
and uvw2 bands (having central wavelengths in units 
of A at 3465, 4392, 2600, 192 8 respectively) with the 
UV-Optical Telescope (UVOT; iRoming et ail [20051 ) on 
board the Swift satellite (Table Q]) . We used the stan- 
dard pipeline reduced image products, co-added and ex- 
posure corrected within the XLMAGE0 environment. For 
the photometry we used the standard 5" aperture. The 
reddening cor rection of E(B — V) = 0.1 77 was applied 
accord ing to iSchlegel et al.l (|1998|) and iCardelli et al.l 
(fl989h . 

Palomar. <?-band (central wavelength 5240 A) obser- 
vations were c arried out with th e robotic Palomar 60- 
inch telescope (|Cenko et al.l I2006T ) (see Table [1]). Data 
were reduced with standard IRAF routines. Photomet- 
ric calibration was performed relative to the USNO-B1 
catalogQ which leads to a systematic contribution to the 
photometric uncertainties of ~ 0.5 mag. 

2.2. X-ray observations 

INTEGRAL. The INTEGRAL satellite observed the 
source during the North Ecliptic Pole Key P rogram with 
its imager IBIS/ISGRI (|Lebrun et al.ll2003fl which oper- 
ates in the 17 - 1000 keV range (Table [J). The observa- 
tions were performed with a rectangular 5x5 dithering 
patter n. We used the standard Off-line Science Analysis 
(OSA; iCourvoisier et al.l I2003D software version 7.0 for 
the ISGRI analysis. Most recent matrices available for 
standard software (isgri_arf_rsp_0025.fits) were used for 
spectral analysis. We used 83 science windows for a total 
amount of 661 ks of exposure time. Data screening was 
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performed according to the median count rate with re- 
spect to each science window and its distribution. After 
data cleaning the effective exposure on the source is 466 
ks. INTEGRAL/SPl upper limits were obtained with 
the SPIMODFIT software, which is available in OSA. It 
performs spectral model fitting for point sources and dif- 
fuse emission b ased on the maxim um likelihood method 
as described in iPetrv et al.l (|2009D . Due to the adopted 
dithering pattern, 1ES 1959+650 is outside the field of 
view of INTE GRAL'S Join t European Monitor for X rays 
(JEM-X: lLu"nd et al.l[200l for at least 77 % of the obser- 
vation time. For the same reason the O ptical Monitoring 
Camera (OMC: [Mas- Hesse et all 120031) could not point 
to the source. Therefore, no flux measurement could be 
obtained from these instruments. 

Swift/XRT The X-Ray Telescope (XRT; 
iBurrows et al.1 120051 ) on board Swift observed the 
blazar in 2007 November (Start date: UTC 2007 
November 30 11:57:33) for an exposure time of 1283 s. 
Data processing, screening and filtering were done using 
the FTOOL xrtpipeline included in the HEAsoft 6.3 
distribution. 

3. RESULTS 

3.1. Variability analysis 

Blazars are strongly variable objects on timesca les as 
short as minutes to hours (e.g.. iFossati et aT1 l2008). The 
long exposures used with IBIS/ISGRI may therefore be 
affected by variations in flux and spectrum. Thus, we 
analyzed the available multiwavelength data set for indi- 
cations of variability. 

The non-simultaneous (to IBIS/ISGRI, Swift/XRT 
and Swift/VVOT; see Table [J) Palomar 60- inch obser- 
vations show a very low variability in brightness of the 
source of 10 % at a significance of ~ 3er. 

The UVOT observations show marginal brightness 
variations by about 10 % on a daily timescale and 20 % 
on inter-day timescales (signific ant at ~ 6a). 

The all-sky monitor (ASM; iLevine et aTlli996D on the 
RXTE) provides regular monitoring of 1ES 1959+650 in 
the 1.5 - 10 keV bandQ. The dwell-by-dwell 90s expo- 
sure light curve does not show flaring activity during the 

http: / /xte. mit.edu/asmlc/ASM. html 
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TABLE 2 

X-RAY SPECTRAL FIT RESULTS. 



Inst. 
(1) 



Fit Mo 

(2) 



(3) 



XRT 
XRT 
ISGRI 
joint 
joint 



ab pi 
ab pi 

Pi 
ab bkn 
log-par 



TToT" 

0.03 
fixed 



(4) 



o.i7: 



2.1 
2.0 



TTTT" 



b 

(5) 



(6) 



Norm 
(7) 



fixed 
fixed 



1-98:1 
i 3° 1 



7 fi 1 - 
'•"0.5 

6° 1 
u - °Q.i 



7 2 3 - 8 
' - z 0.1 



4 -fTJX" 

o J:! 

o -°0.1 
7 S 2 ' 6 

o oO.l 
o 4 0.1 



X 2 /d°f 
(8) 



Flux a 
(9) 



178/192 
196/193 

0.6/3 
200/198 
323/199 



1 n n0.02 

QQ2 

1 nqO-0 2 
i.ua 02 

1 4O.3 



Explanation of columns: (l)instruments; (2)fit model; (3)column density (units of 10 cm -2 ); (4)spectral index; (5)spectral index after 
break; (6)break energy; (7)normalization (units of 10 -4 ph keV - 1 cm -2 ); (8)statistical fit result; (9)flux units are 10 -8 erg cm -2 s -1 for 
XRT and joint spectra and 10 — 11 erg cm -2 s _1 for IBIS/ISGRI spectrum; "fluxes are computed in the 0.6 — 6.0 keV and 20 - 60 keV 
energy ranges for the XRT and ISGRI respectively. 

1ES 1959+650 with a flux of ~ 1.5 x lO" 11 erg cm- 2 
s -1 in the 20 - 60 keV range. The low Galactic re- 
value (1.0 x 10 21 cm" 2 ; iDickev fc Lockmanlfl990[) does 
not affect the IBIS/ISGRI energy range, and we do not 
see evidence for further absorption or deviations from a 
single power law. 

The best-fit result (x 2 red ~ 0.9) for the Swift/XRT spec- 
trum is given by an absorbed power-law model, with ab- 
sorption parameter free to vary. The derived N#-value 
(l-7f\ x 10 21 cm -2 ) is a factor 1.7 larger (at a signifi- 
cance of ~ 2cr) than the Galactic one. 

We do not see significant variability neither on long 
timescales nor on short timescales in the hard X-ray (20 
- 40 keV) and in the X-ray (1.5 - 10 keV) bands dur- 
ing the multifrequency campaign (see Figure [l}. There- 
fore, we combine the relatively short (simultaneous to 
IBIS/ISGRI) Swift/XRT observation (0.6 - 6 keV) with 
the averaged IBIS/ISGRI spectrum. We model the data 
best with a broken power-law model with the absorption 
parameter fixed to the Galactic value. Furthermore, we 
have applied a log-parabolic model: 




F{E)^K{E/E 1 ) 



-(a+b logiE/Et)) 



(1) 



54431 54432 
MJD 



Fig. 1. — Upper panel: ASM light curve of the source in the 
1.5 - 10 keV energy band. Black diamond data points and red 
rectangles represent the dwell- by-dwell 90 s exposure and the one- 
day average light curves respectively. Any flaring activity can be 
excluded. Lower panel: ISGRI light curve shows a constant flux 
level (revolutions 624 and 625 are co-added due to the low exposure 
of the former revolution) . 

INTEGRAL /IBIS observations (Figure [Q upper panel). 
The one-day average light curve excludes the presence 
of any flare, and it shows flux variability (< 20%), well 
within the errors. 

The 20-40 keV IBIS/ISGRI light curve binned per revo- 
lution is shown in Figure [1]- lower panel (revolutions 624 
and 625 are co-added). The low detection significance 
of the source in each individual revolution (~ 5cr) does 
not allow a proper light-curve fitting. If binned to SCW 
level (time binning of ~1 hr), the light curve does not 
allow to rule out some marginally significant variability 
due to the low signal-to-noise ratio (S/N). However, we 
can exclude any flaring activity in the 20-40 keV band. 

3.2. Swift/XRT and IBIS/ISGRI spectra 

The fit results for the single Swift/XRT and averaged 
IBIS/ISGRI spectrum are shown in Tabled The spec- 
tra were fitted using XSPEC 12 and the latest available 
response matrices for calibration. IBIS/ISGRI detected 



This c urved model was first proposed bv lLandau et al.l 
(1986) t o describe the synchro tron component of BL Lac 
objects. iMassaro et al.l (|2004l ) used it to describe the syn- 
chrotron X-ray component of the spectrum of TeV BL 
Lac objects. The model is explained in physical terms 
by means of radiative cooling processes (via synchrotron 
and inverse- Compton) of the high-energy electron popu- 
lation injected with power-law slope. The log-parabolic 
law is a rather simple analytical formula related to physi- 
cal parameters of the source. It is applied to synchrotron 
broadband spectra since it better describes the spectra 
compared to power laws with exponential cutoff. The 
photon index a is considered at energy E\, and the pa- 
rameter b describes the curvature of the spectrum. The 
fit result is reported in Table [2] 

4. THE SPECTRAL ENERGY DISTRIBUTION 

Iii order to construct a multiwavelength spectrum of 
the source we use the data sampled during the simulta- 
neous multifrequency campaign in 2007. At TeV energies 
we show non simultane ous archive HEGRA IACT data 
jAh aronian et al.l 12003). when the source showed a ma- 
jor outburst (in 2002). Such activity at TeV energies 
was not reported during our multifrequency campaign in 
2007. Therefore these data are not used to constrain our 
SED model. 

Given the lack of variability on short and long 
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TABLE 3 

List of parameters used to construct the theoretical SEDs. 



Model R »j osc B D 7 min 7 max q 7 pm i n 7pmax <2 P L p 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 

leptonic 0~5 10 14 19 1x10 s 6xl0 4 L~85 

lcpto-hadronic 2 5 20 19 8xl0 2 4.5xl0 4 1.9 lxlO 3 1.2xl0 9 1.9 3.5 

Explanation of columns: (1) = SED model; (2) = radius of emitting region in units of 10 14 cm; (3) = escape time parameter: t esc = r\ X R/c; 
(4) = magnetic field in Gauss; (5) = Doppler factor; (6) and (7) = minimum and maximum random Lorentz factors of the injected electrons; 
(8) = slope of the injected electron distribution; (9) and (10) = minimum and maximum random Lorentz factors of the injected protons; 
(11) = slope of the injected proton distribution; (12) = kinetic power of relativistic protons in units of 10 46 erg s — 1 . 




Log v [Hz] 



Fig. 2.— Observed SED of 1ES 1959+650 in 2007 November and 
December. The model fits to the SED (see the text) are overplotted 
(black solid and orange solid curves are the leptonic and lepto- 
hadronic models respectively). 77 absorption by the extragalactic 
background light (EBL) is accounted. 



timescales, we used the average values in each energy 
band to construct the SED. The SED of 1ES 1959+650 
shows the typical two-component structure of an HBL 
(Figure [5]). The synchrotron peak can be constrained 
by the XRT spectrum (a ~ 2, see Table [5]), whereas 
the hard X-ray-soft 7-ray data constrain the onset of 
the high-energy component. Indeed, the log-parabolic fit 
model for the synchrotron emission is not able to repro- 
duce the jointly fitted XRT and IBIS/ISGRI data (see 
fit result in Table [5]). This, in turn, suggests that the 
IBIS/ISGRI spectrum represents the onset of the high- 
energy component. The high-energy peak is very poorly 
constrained due to the lack of simultaneous 7-ray data. 
During its first year of survey, the Fermi Large Area 
Telescope (LAT) d etects the source in quiescence state 
(jAbdo et al.l I2010l ). Therefore the power output mea- 
sured by LAT is not comparable to the power output 
in the same energy range inferred with our model that 
refers to an active state of the source. The LAT measure- 
ment is a factor ^6 lower. The Fermi-L AT observations 
are taken ~1 year after our multifrequency campaign of 
2007. 

The data are modeled by a pure leptonic SSC model 
(see Figure [5] black s olid line) , using the equilibr ium ver- 
sion of the code of iBottcher fc Chiang! (120021) , as de- 
scribed in more detail in (jAcciari et al.H2009f ). The ge- 
ometry of the emitting region is a spherical volume V b ' of 
radius i?b in the comoving frame. It moves with respect 
to the observer with a bulk Lorentz factor T (speed /?rc) 
at an angle o bs, resulting in relativistic beaming deter- 
mined by the Doppler factor. Ultrarelativistic leptons 
are injected into the emission region with a power-law 



distribution (in the comoving frame): 

Q[ ni (i;t) = Q i o\t)"f- q [cm-h- 1 ] for 71 < 7 < 72 

where the normalization is determined by the injection 
power ij n j. The code finds a self-consistent equilibrium 
between particle injection, radiative cooling due to syn- 
chrotron and synchrotron self-Compton (SSC) losses and 
particle escape on a timescale t csc — r/ csc R/c. The ef- 
fect of 77 absorption by the extragalactic background 
light (EBL) is tak en into account using the model of 
iFinke et~al . (2010). The fit is constrained by the syn- 
chrotron component from the optical to X-rays as well 
as the onset of the SSC component at hard X-rays - soft 
7-rays. The hard synchrotron spectrum as well as the 
unusually high level of the low-frequency end of the SSC 
component (see Figure [5]) requires the choice of a very 
hard injection spectrum with a slope of q — 1.85. The 
optical data point from the Palomar observations were 
taken up to 10 days after the X-ray and UV pointings 
with XRT and UVOT respectively. They cannot be rec- 
onciled with our model SED. This is not surprising, as 
variability by a factor of ~ 2 on timescales of w eeks is 
not uncommon in this object ([Villata et al J 12000). The 
parameters used for the fit shown in Figure [2] are listed 
in Table [3J The equilibrium particle distribution found 
by the code corresponds to a kinetic power in relativistic 
electrons of L e = 8.5 x 10 42 erg s _1 , while the mag- 
netic field of B = 14 G yields to a power in Poynting 
flux of Lb — 7.4 x 10 41 erg s _1 . Hence, the magnetic 
field energy density is a factor es = Ls/L e = 0.09 
below equipartition. The hard injection index of q = 
1.85 (constrained by the mere synchrotron component) 
is inconsistent with standard first-order Fermi acceler- 
ation at relati vistic shocks, which predicts an index of 
q ~ 2 .2 - 2.3 (jAchterberg et al.ll200ll lEllison fc Double! 
I2004f h This might indicate a substantial contribution 
to par ticle acceleration from second-order Fermi acc eler- 
ation (IVirtanen fc Vainiol [20051 : IStecker et al.l 127)01 . A 
similar conclusion was also reached when modeling the 
very hard X-ray and Fermi 7-ray spectrum of the HBL 
RGB J0710+591, recently detected at very- high-energy 
(VHE) 7-rays by VERITAS. 

As an alternative to the pure leptonic model, we have 
applied a semi-analytical lepto-hadronic model shown by 
the orange solid line in Figure [U This model assumes, 
in addition to a leptonic component similar to the one 
used for the leptonic model described above, a power- 
law population of relativistic protons extending out to 
energies beyond the threshold for P7 pion production 
on the electron-synchrotron radiation field. The produc- 
tion rates of final decay products (electrons, positrons, 
7T° decay photons, and neutrinos) are calculated using 
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the analytical templates of Kelner fc Aharoniartl (|2008T) . 
Synchrotron emission of secondaries is calculated using 
a 3u{l) K i/ 1 / 3 e~ !y / l/ °( 7 ) approximation. The ir° decay 
photons as well as synchrotron emission from the first- 
generation pairs from charged pion decay are produced 
predominantly at ^> TeV energies, at which the emission 
region is highly opaque to 77 absorption. Therefore, the 
radiative power at those energies is redistributed to lower 
frequencies through electromagnetic cascades. We em- 
ploy a semi-an alytical tr e atme nt of the cascading process 
as described in Bottcher (2010). The inferred parameters 
are reported in Tabel El and they are in good agreement 
with the leptonic model. 

5. CONCLUSIONS 

We have performed a multiwavelength campaign on 
1ES 1959+650 in 2007. This is the first and only de- 
tection by IBIS/ISGRI of the source. The source was 
monitored nearly simultaneously from optical to hard X- 
ray energies. In the optical to UV band the most evi- 
dent flux variations are of the order of 20 % on inter-day 
timescales. At hard X-rays the source flux stayed con- 
stant. The compiled SED allowed us to derive the phys- 
ical parameters of the source. The SED clearly shows 
the simultaneous detection of both the synchrotron and 
the rise of the high-energy (inverse-Compton) emission 
components. This is the first time that the transition 
region between the synchrotron and inverse-Compton 
components is detected for 1ES 1959+650, and has so 
far been measured only for a few other BL Lac ob- 
j ects of intermediate ty pe. The best example is ON 231 
(Tagliafcrri ct al. 2000) which is an Intermediate BL Lac 
object (IBL). These objects are characterized by SEDs 
peaking at frequencies intermediate to LBL and HBL. 
Other examples , not as clear as the previous one, ar e 
PKS 2155-304 (IKubo et alJU99l IFoschini et all [ 2008). 
another IBL and S5 0716+ 714 (IFoschini et all 120061 : 
iGiommi et aJ]|1999T ). an LBL. IFoschini et all (|2008l ) sug- 
gest that the source is in a continuous high active state 
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